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Abstract 

Solar neutrino oscillations with wavelengths comparable to the Earth-Sun dis¬ 
tance provide a viable explanations of the long-standing solar neutrino deficit. 
They imply a time-dependent modulation of the solar neutrino flux due to 
the eccentricity of the Earth orbit. Motivated by this testable prediction, we 
propose a Fourier analysis of the signal observable in real-time, solar neutrino 
experiments. We give the general expressions of the Fourier coefficients and 
of their correlated uncertainties in the presence of background. The expres¬ 
sions assume a particularly compact form in the case of two-flavor neutrino 
oscillations in vacuum. We discuss the sensitivity to the lowest harmonics of 
the new-generation, high-statistics experiments SuperKamiokande, Sudbury 
Neutrino Observatory, and Borexino. 
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The phenomenon of neutrino oscillations in vacuum [|l| represents a possible solution 
to the long-standing solar neutrino problem Considering for simplicity only the 
hrst two neutrino families v^. and the v^. survival probability P at a distance L from the 
Sun is given by 


P(E) 



— COS 


‘IE ) 


( 1 ) 


where E is the neutrino energy, 6m? is the neutrino squared mass difference, and 9 is the 
vacuum mixing angle. 

A tentative evidence for the oscillating term in Eq. (|l|) comes from the P-dependence 
of the solar neutrino dehcit [Q as inferred from the four pioneering solar u experiments [0. 
The ellipticity of the Earth’s orbit implies a further, striking signature of the oscillation 
phenomenon, namely, the L-dependence of the observed flux (in addition to the trivial 1/L^ 
geometric factor) pj. In particular, the survival probability P in Eq. (|^ is modulated in 
time by the periodic variation of L which, at hrst order in the eccentricity {e = 0.0167), is 
given by 


L{t) = To ^ cos -h 0(e^) 


( 2 ) 


with Lq = 1 AU, T = 1 yr, and t = 0 at the perihelion. The hrst-generation experiments 
have not collected enough statistics to test the L-dependence of the solar z/ hux |^. New- 
generation experiments should instead be able to probe the structure of the neutrino signal 
in the time domain ra¬ 


in this work we propose to study the neutrino signal in the frequency domain through 
a Eourier analysis of the periodic variations associated to havor oscillations in vacuum. 
This approach is particularly suited with real-time, high statistics experiments such as Su¬ 
per Kamiokande (operating), the Sudbury Neutrino Observatory [Q (SNO, in construc¬ 

tion), and Borexino (in construction). It will be seen that the different sensitivity of each 
experiment to the lowest harmonics can help in discriminating the value of 6m‘^ (should the 
vacuum oscillation solution be confirmed). Compact expressions for the Fourier coefficients 
and for their uncertainties will be given. 
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In general, the neutrino event rate R{t) at the time t is the sum of a signal S{t) and of 


a (supposedly constant) background B, 


R{t) =B + Sit) . 


( 3 ) 


For symmetry reasons, the analysis can be restricted to the time interval [0, T/2]. It is 
understood that events collected in subsequent half-years must be symmetrically folded 
in this interval. The data sample consists then of N events collected at different times 
{ti}i<i<N, with ti G [0, T/2] and N equal to the total sum of background and signal events, 
N = Ns + Nb- Notice that, in general, the background-signal separation can be performed 
on the average rates, but not on an event-by-event basis. 

The expansion of the signal in terms of Fourier components fn is dehned as 

SW = s('l + 2f;/„cos“) , (4) 

where S is the time-averaged signal 


2 NR 

S = -j^ dtSit) 


( 5 ) 


The n-th harmonic corresponds to a period of 1/n yr. The explicit form of fn reads 


rT/2 


fn = ^ dt R{t) COS 
bl Jo 

1 27rnti 

= -TT- / COS ——— 

Ns h T 


2Tint 
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( 6 ) 

( 7 ) 


where Eqs. (|]) and ( 0 ) represent the theoretical dehnition and the experimental determina¬ 
tion of the /n’s, respectively. Although in Eq. (|) one can replace i?(t) with S'(t), the use of 
Rit) is more general, since background events do contribute to the statistical uncertainties. 
Notice that the sum in Eq. (0) does not require an event-by-event separation of signal and 
background, and does not involve any binning in time. 

Each fn is a linear combination of Poisson random variables = R{t)dt [Eq. (|^)], which 
represent the total number of events collected in the interval [t, t + dt]. If only statistical 
errors are considered, then var(^t) = and cov(^t, N) = 0, since fluctuations at different 
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times t and t' are uncorrelated |T^. The (co)variances of linear combinations of independent 
random variables are given by []^ var(X]t var(,^t) and covij^t A'6') = 

atAvar(,^i). It follows that 

r'^/2 ^ / 2 27mf\ ^ 


/ . N r/% T./ N / 2 27mtY 

var(/„) = dt R{t) i^— cos , 


rT/2 

COv{fn, fm) = / dtR{t) 

Jo 


4 27mt 27rmf 
cos —r— cos ■ 


(5T)2 T T 

The hnal result for the one-sigma error cr„ affecting /„ and for the correlation p^n ^ n) 

is: 


^ n. 


1 + /2n + B / S 


ST 


( 8 ) 


Pmn — 


fm+n “1“ /*|m—n| 


( 9 ) 


(1 + / 2 m + B/ S')(l -|- /2n + -B/ S) 

The values of (T„ and pmn can be expressed in terms of measured quantities with the substi¬ 
tutions ST/2 —i> Ns and B/S ^ Nb/Ns- 

In the standard (std) case, i.e. in the absence of oscillations, the Fourier transform of the 
signal 

/ r \ 2 


is trivial, only the hrst coefficient being nonzero and equal to the eccentricity £, 


fn^ = ednl . 


( 10 ) 


( 11 ) 


Moreover, in the standard case the statistical errors do not depend on n. 


std ^ l ^B + Ns 

V 2iV| 


( 12 ) 


and thus form a “white noise” affecting all the harmonics. The error correlations, as derived 
from Eq. (|), are given by = {eNs/N)5m,n±i < ^ and thus are practically negligible. 

In the case of two-flavor oscillations the signal is proportional to 


4 




















(13) 


Sit) IdE ME) \aME)P + ff„(B)(l - P)|, 

where \{E) is the neutrino energy spectrum, Ue{E) and <J^{E) are the z/g and interaction 
cross sections, and the probability P is given by Eq. (|l|) with L as in Eq. (H). (It is understood 
that (Te and cr^ are corrected for energy threshold and resolution effects in each detector.) 

Given the signal S{t) as in Eq. ([I^), the time integration in the expression of the Fourier 
coefficients 


fn = 


S (t) cos 


2nnt 

T 


iPdtS(t) 

can be performed analytically. The hnal result can be cast in the following, compact form: 

eSni — sin^ 26* Dn{5m?) 


fn = 


1 — sin^ 26 

where the detector-dependent functions are given by 


(^ > 1 ) , 


(14) 


Dn{6m^) = 


2 _ JdEX{ae - cr^) 


(n > 0) 


(16) 


2 / dE X (Tg 

and the universal (i.e., detector-independent) functions Un are given by 

Un{z) = COS (^Z - Jn{ez) , 

where = 5m?LQ/2E and is the Bessel function of order n |^. Notice that, although 
our calculations are of 0{e), all orders in are kept, since ^ may be large. 

We apply now the Fourier analysis to the signal expected in the SuperKamiokande, SNO, 
and Borexino experiments. The SuperKamiokande and Borexino cross sections {u + e ^ 
u' + e) are taken from |]^, and the SNO cross section (z/g -t- d p + p + e) from |1l8[ . 
These cross sections are corrected for energy threshold and resolution effects (see, e.g., [0). 
In particular, we consider a prospected threshold of 5 MeV for the recoil electron kinetic 
energy Tg in SuperKamiokande and SNO, and an analysis window Tg G [0.25, 0.8] MeV for 
Borexino. The energy resolution is assumed to be 16% at 10 MeV for SuperKamiokande, 
11% at 10 MeV for SNO, and 8% at 0.5 MeV for Borexino, scaling as for different recoil 
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energies. The Fourier components of the vacuum oscillation signal are calculated through 
Eq. ([I^) in the 6m‘^ range relevant for the solution of the solar neutrino problem and then 
plotted in Fig. 1. 

Figure 1 shows the deviations of the Fourier components from their standard values, 
/„—as functions of for maximal mixing (sin^ 29 = 1). The deviations decrease with 
decreasing mixing (not shown), although not exactly in proportion to sin^ 29 [see Eq. (0)]. 
The three subfigures refer to SuperKamiokande, SNO, and Borexino (top to bottom). The 
gray, horizontal bands correspond to the Tier statistical error relative to the standard (no 
oscillation) case for Ns = 10^ events and Nb = 0 (cr„ = 0.0071). Error bands for different 
choices of Ns or Nb can be obtained through Eq. (|^). 

It can be seen from Fig. 1 that the SuperKamiokande and SNO experiments are sensitive 
only to the first harmonic, i.e. to variations of the signal with a period of 1 yr. The amplitude 
of the second harmonic is, in fact, much smaller than the ±lcr error in both cases. This is in 


agreement with the results in |^, where the SNO and SuperKamiokande signals were shown 
to vary almost sinusoidally in time. A single parameter (/i) is then sufficient to characterize 
the results of each of these two experiments. For 10^ events and no background. Fig. 1 shows 
that the genuine n = 1 oscillation component (fi—e) may reach a significance of ~ Scr ( 7a) 

in the most favorable case for SuperKamiokande (SNO). For a signal-to-background ratio 
equal to one {Nb = Ns), the statistical signihcance of the vacuum oscillation signal would 
be lower by a factor \/2. This implies that the detection of an unmistakable signal for 
vacuum oscillation in the time or frequency domain will require very high statistics and 
good background rejection in both experiments. Even in the best conditions, there are some 
ranges of 6m‘^ where the expected time-modulation of the signal is small and undetectable 
(e.g., for 6m? ~ 1.6 x 10“^° eV^ in Fig. 1) |1imp0| . In such ranges, however, independent 
oscillation signals might show up in the energy domain as distortions of the recoil electron 
spectrum 00- 

The situation is much more favorable for the Borexino experiment, since its signal is 
dominated by the monoenergetic ^Be solar neutrinos {E = 0.862 MeV), while the Su- 
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perKamiokande and SNO signals are smeared by the broad ®B neutrino spectrum {E ^ 15 
MeV). The third panel of Fig. 1 shows, in fact, that at least one of the first three Fourier 
coefficients is larger than the (representative) statistical error in the 5m? range of interest. 
Moreover, when one of the harmonics is small, another is large (and vice versa), so that 
no “holes” are left in the sensitivity to the 5m? variable, provided that (at least) two of 
the first three Fourier components are measured. Notice also that the relative amplitude of 
the Fourier coefficients in Borexino is strongly dependent on 5m?, and thus the detection of 
two nonzero harmonics would be of great help in discriminating a preferred range of 5m?. 
The relative amplitude of the first Fourier coefficient in SuperKamiokande (or SNO) and 
Borexino is also dependent on 5m? and, therefore, the combination of all the experiments 
will enhance the resolution in 5m?. 

In conclusion, we have performed a Fourier analysis of the signal expected in the Su- 
perKamiokande, SNO, and Borexino solar neutrino experiments. This method fully exploits 
the real-time features of the three detectors and requires no binning in time. Expressions for 
the correlated uncertainties of the Fourier components have been worked out in the general 
case and, in particular, for no oscillation. Compact expressions for the Fourier coefficients in 
the presence of 2v oscillations have been given. The method has been applied to the analy¬ 
sis of the signals expected in the new-generation experiments. SuperKamiokande and SNO 
are shown to be sensitive only to the hrst harmonic, while Borexino is also sensitive to the 
second and third harmonic. The combination of different Fourier component measurements 
is highly selective in 5m?. 

We thank G. Bellini, M. G. Giammarchi, and A. lanni for useful correspondence, and P. 
I. Krastev for helpful suggestions. 
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FIGURES 


FIG. 1. Deviations of the Fourier coefficients fn from their standard values = e5ni, as func¬ 
tions of dm? for maximal 2u mixing (sin^ 29 = 1). In each of the three panels (SuperKamiokande, 
SNO, and Borexino), the gray, horizontal band represents the ±la statistical uncertainty associ¬ 
ated to the standard predictions for 10^ events with no background. SuperKamiokande and 
SNO appear to be sensitive only to the first harmonic (solid curve), while Borexino is sensitive also 
to the second and third harmonic (dashed and dotted curve, respectively). Notice that the vertical 
scales are different, but the absolute width of the gray error band (an = ±0.0071) is the same for 
the three experiments. 
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FIG. 1 Deviations of the Fourier coefficients /„ from their standard values = £6ni, as func¬ 
tions of 6m? for maximal 2iy mixing (sin^ 29 = 1). In each of the three panels (SuperKamiokande, 
SNO, and Borexino), the gray, horizontal band represents the ±1(T statistical uncertainty associ¬ 
ated to the standard predictions for 10^ events with no background. SuperKamiokande and 
SNO appear to be sensitive only to the first harmonic (solid curve), while Borexino is sensitive also 
to the second and third harmonic (dashed and dotted curve, respectively). Notice that the vertical 
scales are different, but the absolute width of the gray error band (cr^ = ±0.0071) is the same for 
the three experiments. 
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